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Abstract. In this work the analysis of mixing parameters of the system involving 
77, rl mesons and some third massive state G is carried out. We use the generalized 
mixing scheme with three angles. The framework of the dispersive approach to 
Abelian axial anomaly of isoscalar non-singlet current and the analysis of exper- 
imental data of charmonium radiative decays ratio allow us to get a number of 
quite strict constraints for the mixing parameters. The analysis shows that the 
equal values of axial current coupling constants /g and /o are preferable which 
may be considered as a manifestation of SU (3) and chiral symmetry. 



1 Introduction 

This work is developing the approach of the papers 13115^ and is devoted to the significant 
problem of mixing of pseudoscalar mesons. It is especially important with a number of current 
and planned experiments. 

The problem of 77-77' mixing has been studied for many years. The usual approach with one 
mixing angle dominated for decades, but in the recent years the more elaborated schemes appear 
to be unavoidable [SHS]. In particular, the theoretical ground of this was based on the recent 
progress in the ChPT [5HTT]. On the other hand, it was shown, that the current experimental 
data cannot satisfactory describe the whole set of experiments within the one-angle mixing 
scheme. 

The mixing schemes are usually enunciated either in terms of SU (3) or quark basis. In 
our paper [5] we construct and use the generalization of 5J7(3) basis similar to the mixing of 
massive neutrinos. This is because we use the dispersive approach to axial anomaly ( |12j . [13] 
for a review) to find some model-independent and precise restriction on the mixing parameters. 

It was shown that any scheme with more than one angle unavoidably demands an additional 
admixture of higher mass state. If we restrict ourselves to only one additional state G (denoted 
as a glueball without really specifying its nature) then the general mixing scheme can be 
described in terms of 3 angles. In particular cases the number of angles can be reduced to two. 

In the paper [2] the analysis of different conventional (and most physically interesting) 
particular cases was performed (including two-angle mixing schemes) basing on the dispersive 
representation of axial anomaly from one side and charmonium decays ratio from the other 
side. 
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The main conclusion of the paper [5] is that in aU considered cases the only reasonable 
solutions appear at fg, — fo — /tt- The main aim of this work is to check whether this relation 
remains valid in the most general case with some specific constraints imposed. 

This paper is organized as follows. In the Sec. [5] we introduce our notation and the general 
approach to the mixing. In Sec. [3] we derive the basic equations relying on the dispersive 
approach to Abelian axial anomaly of isoscalar non-singlet current and the charmonium 
radiative decay ratio Rj/^, while in Sec.0]we perform the numerical analysis of these equations. 
Finally, in Sec. Owe present the conclusion. 



2 Mixing scheme 



We start with a (A^-component) vector of physical pseudoscalar fields consisting of the fields of 
the lightest pseudoscalar mesons and other fields: 



G 



(1) 



We are not able to specify the physical nature of the other components with higher masses, 
the lowest of which G can be either a glueball or some excited state Q- Let us also introduce. 



following [IIKIT], a set of SU{3) fields ip3, 
(sterile) fields gi {<Pi,i = 4. .TV) 



, ipo <?2, ^^s) and complement them with other 



.9 



(2) 



The three upper fields ip3,(ps, Vo are the only ones which define the generalized PCAC 
relation for axial current J"^ = 97^^75 (no summation over a contrary to j and k is assumed): 



^^. JU = = FajM.k'Pk, a - 3, 8, 0, j, k = l..iV, 



(3) 



where AC is the mass term in the effective Lagrangian with a non-diagonal mass matrix M (as 
fields <Pk are not orthogonal to each other): 



AC = -*'^M*, 
2 



and F is a matrix of decay constants Q: 



\f3 
F = I /8 
/o 



(4) 



(5) 



In order to proceed from initial SU{3) fields # to physical mass fields $ the unitary (real, 
as the CP-violating effects are negligible) matrix U is introduced 



^ Note, that the mixing with the excited states is usually(e.g. |14II15) 1 supposed to be suppressed. 

^ Note, that matrix of decay constants F is non-square expressing the fact that generally the number 
of SU{3) currents is less then the number of all possible states involved in mixing. The similar situation 
takes place (see e.g. [18]) in one of the extensions of the Standard Model - neutrino mixing scenario 
involving sterile neutrinos. 
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* = U* (6) 

that diagonalizes the mass matrix 

UMU'^ = M = diag{mlo,m^^,m^,,ml;, . . .), (7) 

where 771^^, m^, mj^i and iuq are the masses of the tt, 77, 77' mesons and glueball state G, 
respectively. 

Simple transformations of Eq.(j3]) read: 

9^J^5 = FU^Md (8) 

This formula is close to those obtained in [TBUTT] (in the limit of small mixing). When the decay 
constants are equal, it is reduced to formula (3.40) in [TH]. 

The matrix elements of S^J^s between vacuum state and physical states \'Pk) 

{0\d^j;^\¥k) = FnU^Mfk (9) 
can be compared to the standard definition of the "physical" coupling constants of axial cur- 
rents* 

(0| J^sl^) = ^/fe^P- (10) 
From © and ([TUI) follows the relation 

f^^F^{U^)l^UU^)l. (11) 

This expression (recall, that there is no summation over a ) clearly shows that are obtained 
by multiplication of each line of U"^ by respective coupling and form a non-diagonal (contrary 
to F) matrix. 

Taking into account the well-known smallness of tt*' mixing with the r/, ry' sector |16 1 I17| , [20] 
and neglecting all higher contributions we restrict our consideration to three physical states 
•q, rj', G and two currents J^g, J^g. Then the divergencies of the axial currents (recall, that G is 
a first mass state heavier than 77'): 



9,.JU\ _ (h 




(12) 



Exploring the mentioned similarity of the meson and lepton mixing, we use the Euler 
parametrization for the mixing matrix U (we use notation Ci = cosOi, Si = sinOi): 



(C8C3 - C0S3S8 -C3S8 - C8C0S3 S3S0 
S3C8 + C3C0S8 -S3S8 -I- C3C8C0 -C3S0 I . (13) 
ssso CsSo Co 

In the following consideration we will need the divergency of the octet current d^J^^, so let 
us write it out explicitly: 

9tJ.J^5 = fsimlvicscs - CqSsSs) + m^yr/{s3Cs + c^coSs) + m^GissSo)). (14) 

As soon as in the chiral limit J^g should be conserved, from Ea. lfTHl follows that coefficients 
of the terms m^,,?7ig. must decrease at least as {nij^ / rrirj' ,g)'^ ■ More specifically, we expect the 
following limits for the terms of Eq. lfT^ : 

M <(!I^)\ (15) 



IS3C8 + C3C0S8I V™G 
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3 Abelian axial anomaly and charmonium decays ratio 

In our paper the dispersive form of the anomaly sum rule will be extensively used, so we remind 
briefly the main points of this approach (see e.g. review [T3] for details). 

Consider a matrix element of a transition of the axial current to two photons with momenta 
p and p' 

T^a^{p,p') = (p,pVm5|0) . (16) 

The general form of T^ap for a case p^ = p'^ can be represented in terms of structure 
functions (form factors): 



T^afsiP,?') = Fi{q'^)qf,€al3paPpP'„ + 

1 TP ( 2^,Pa , P'l3 

7.F2{q )[-^^nPp<yPpPa- - 
z, p p 



F2{q'^)[^epi3pcrPpP'a- " -^^^.i^apaPpPa " E/iQ/Ja (p ^ p')<t] , (17) 



where q = p + p' . The functions Fi{q'^), F2{q^) can be described by dispersion relations with 
no subtractions and anomaly condition in QCD results in the sum rule: 

oo 

j Im Fi{q^)dq^ = 2aNcJ24 ^ (1^) 



where Cq are quark electric charges and Nc is the number of colors. This sum rule [21j was 
developed by Jiff Hofejsi [55], and later generalized [33]. Notice that in QCD this equation 
does not have any perturbative corrections [24], and it is expected that it does not have any 
non-perturbative corrections as well due to the 't Hooft's consistency principle [5S]. It will be 
important for us that as ^ cx) the function ImFi{q^) decreases as l/(?^ (see discussion 
in Ref. ^2^). Note also that the relation (fT8|) contains only mass- independent terms, which is 
especially important for the 8th component of the axial current J^g containing strange quarks: 



1 



Jt,5 = -^iuip.l5U + d-f^-y^d - 2s7^75s) . (19) 
The general sum rule (|T8| takes the form: 

oo 

j ImF,{q')dq^^^a{el + el-2el)N, = ^^a, (20) 



where = 2/3, ed — — —1/3, Nc = 3. 

In order to separate the form factor Fi{q^), multiply TpapiPjP') by q^ilq^ ■ Then, taking the 
imaginary part of Fi{q^), using the expression for dp_J^^ from Eq. (jl2p and unitarity we get: 

ImF^{q')=Imq^^{2j\jll^\0)^ 

- ^(27 I [m'f^i]{csC3 - cosgsg) + mfy'q'isscs + cgCQSg) + m'^^GssSo] \ 0) = 

'^f8[Ar,S{q'^ - rn^)(c8C3 - Cqs^ss) + A^,6{q'^ - m^,)(s3C8 + csCoSg) + AoSiq'^ - m'^){ssSo)]. 

(21) 



If we employ the sum rule (|20|) , we obtain a simple equation: 

(cgca - CqS^Ss) + Pisscs + CgCoSg) + 7(s8So) = (22) 

where 
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Ag 



3 ' 





1 




27 /I 



r;— >27 



647r " 



(23) 
(24) 



Note that if we include higher resonances in this equation, they will be suppressed as 1/ m^es 
by virtue of the mentioned above asymptotic behavior of Fi{q^) (x For the last two terms 

in (22) we can specify this constraint as follows: 



|S3C8 + CsCqSsI ^ l\mG ) 



(25) 



As an additional experimental constraint we use, following [26J[27], the data of the decay 
ratio Rj/^ = {r{J/^) ^ 77'7)/(r(J/^') ^ 7^7). 

As it was pointed out in [5S], the radiative decays J/'F — > v{v')l ^-re dominated by non- 
perturbative gluonic matrix elements, and the ratio of the decay rates Rj/xp — {r{J/''P) — > 
i]'j)/{r{J/'I') — > 777) can be expressed as follows: 



Rj/^ — 



{0\GG\ i) 



(0 I GG I ^) 



Pv 



(26) 



where Priin') — — ■m'^f^^y^/Mj^^)/2. The advantage of this ratio is expected smallness of 

perturbative and non-perturbative corrections. 

The divergencies of singlet and octet components of the axial current in terms of quark 
fields can be written as: 



9^^^5 = ^(m„U75U + nidd'-fsd - 2msS'-f5s), 



(27) 
(28) 



Following |26j . neglect the contribution of u- and d- quark masses, then the matrix elements 
of the anomaly term between the vacuum and 77, 77' states are: 



^(0 I GG I 77) ^ (0 I d,J%^ I v) + ^(0 I d.jlf I 77), 
^(0 I GG I 7/) = (0 I d,/;,^ \rj') + ^{0\ d.jjS I V')- 



Using Eq. ((ni), dig), (HH), dSni) we deduce: 



(29) 
(30) 



-Rj/f — 



/o(-s3S8 + C3C8C0) + yj/slsacs + cacoss) 
/o(-C3S8 - C8C0S3) + -^fsicsca - cqSsSs) 



4 / \ 3 



Pv 



(31) 



4 Analysis 



For further analysis it is convenient to rewrite the equations (P2|) . (PT|) in terms of angles 
01 = 98 + 93, 0s and 9o: 



^(Ci + C2 - Co(c2 - Ci)) + ^(Si - S2 + Co(si + S2)) + 7(sgSo) = 



(32) 
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R 



J/1' 



/o(ci - C2 + Co(ci + C2)) + ^/8(si - S2 + co(si + S2)) 
/o(-Sl - S2 - Co(si - S2)) + ;^/8(ci + C2 ~ Co(c2 ~ Ci)) 



Pri' \ 



, (33) 



where 6*2 = 26*8 - Oi. 

The angles ^i, ^g, Oq have the explicit physical meaning. From the definition of the 
mixing matrix U one can see that the angle 9i describes the overlap in the rj — rj' system with 
an accuracy ^ 6q/2 and coincides with their mixing angle as 6*0 0. At the same time Oq is 
responsible for the glueball admixture to — 77' system, and ssSq describes the contribution of 
the glueball state G to the octet component of axial current dJf^^ only. 

In the further analysis we will use the following assumptions: 

I) As we discussed in Sec. [3 the last term in (fT4)) should be suppressed as (TO^/mc)^. So 
we impose the following constraint: 



|g8go| < / rn^ 
\S3C8 + CsCqSsI V™G 



(34) 



II) In sec 3 we found another constraint, which follows from the asymptotic behavior of 
ImFi (see[ 



|s8So| ^ P f rrin 



\S3C8 + C3C0S8I 1 \mG 



(35) 



III) In our numerical analysis we suppose that 7 cannot exceed 1 (i.e. rQ^2-y/'niG ^ 
r,j^2i / Tnfj) . This restriction corresponds to the assumption that 2-photon decay widths of 
pseudoscalar mesons grow like the third power of their masses, or in other words, the glueball 
coupling to quarks is of the same order as for the meson octet states. 

IV) We accept that the decay constants obey the relation /8 ^ /o ^ /tt (for various kinds 
of justification see, e.g., [Slin])- 

For the purposes of numerical analysis, the values of Rj/^ [Rj/i — 4.8 ± 0.6), masses 
and two-photon decay widths of rj, r]' mesons are taken from PDG [22] • Using the values 
niri , rriri' , r^-^2-y , r,y , we see that the relation for the constraint (|34p is more strict than the 
constraint ([55)) . Supposing the minimal mass of the glueball to be of order ma — 3m,; ~ 1.5 
GeV, we get the estimation: 

|S8S0|/|S3C8 + C3C0S8I < 0.1. (36) 

On Fig. [1] the plots of the equations ([32]) and (pS)) in the parameter space {9s, Oi) are shown 
for different values of decay constants fg, /o and mixing angle 9q. The dashed curves denote 
experimental uncertainties. The intersection points of the curves represent the solutions of both 
equations (f32j) . ([33 ]) . The filled area indicates the region, where the constraint ([36]) is valid. The 
plotted range of angle 9i is limited to the physically interesting region, where the solution for 
relatively small angles 6*0 exists. Let us note for completeness, that there is another solution for 
9i ~ 90°, 9o > 50° which does not seem to have a physical sense. 

The numerical analysis shows, that the solution of the equations ([5^ and ([55]) satisfying 
the mentioned above constraint is possible only for rather small mixing angle 9o and for decay 
constants /s, /o close to each other and cl ose to fyri f or /s/Ztt = /o/Ztt = 1.0 the possible range of 
mixing angle 9o is 9o = (0-^25)° (see Fig. |l(a)] 1(c) for demonstration), for fs/frr = JqI j-K = 1-1 



the possible range of mixing angle 9q is = (0 20)° 

There is no solutions for decay constant values fs/ f-n = 1-1, /q/Ztt = 1.0 for any 9q (see Fig. 
l(d)||l(f)" for demonstration), and for any /o < fs in case of f^/ f-n > 1-2. The obtained results 



are quite stable: even if we relax the constraint p6p making its r.h.s. several times larger, all 
the conclusions are preserved. 

Note finally, that this result is in contradiction with the prediction for the decay constant 
/s//^ === 1.34 [TO] obtained in the Large ChPT. 



Will be inserted by the editor 



7 




-150-100 -50 50 100 150 -150-100 -50 50 100 150 -150-100 -50 50 100 150 

(^a [degrees] fe [degrees] fls [degrees] 




-150-100 -50 50 100 150 -150-100 -50 50 100 150 -150-100-50 50 100 150 

fla [degrees] % [degrees] 6a [degrees] 



(d) (/8,/o) = (1.1, 1.0)/. A = (e) (/8,/o) = (1.1, 1.0)/. ,^0 = (f) (/8,/o) = (l.l,1.0)/.,eo = 
0° 5° 30° 

Fig. 1. The solutions of the Eq. p2p (thin curves, blue online) and (|33|) (thick curves, red online) with 
the experimental uncertainties (dashed curves) for different values of the parameters /g, fo and 9o. The 
shaded area indicates the region, where the relation p6p is valid. 

5 Conclusion 

In this paper we studied what can be learnt about the mixing in the pseudoscalar sector from 
the dispersive approach to axial anomaly. 

Our analysis shows that the equal values of axial current coupling constants /§ and /o are 
favorable which may be considered as a manifestation of SU(3) and chiral symmetry. Moreover, 
with a less definiteness the relation ~ /s ~ /o [2] is also supported. 

The analysis demands /§ < 1-2 f.^ which deviates at 10% level from the results of calculations 
within the chiral perturbation theory (/§ — 1.34/7r) |10j . 

The value of the mixing angle 6q, which is responsible for the glueball admixture to the 77— 
is limited to 00 < 25° for (/8,/o) = (1.0, 1.0)/^ and to 60 < 20° for (/8,/o) = (1.0, 1.0)/^. 

The improvement of the experimental data of Rj can significantly limit the constraints 
for the parameters 6*0, and fs, fo- 
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